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Abstract We model the gravitational behaviour of a radiating star when the exte- 
rior geometry is the generalised Vaidya spacetime. The interior matter distribution 
is shear-free and undergoing radial heat flow. The exterior energy momentum ten- 
sor is a superposition of a null fluid and a string fluid. An analysis of the junction 
conditions at the stellar surface shows that the pressure at the boundary depends 
on the interior heat flux and the exterior string density. The results for a relativis- 
tic radiating star undergoing nonadiabatic collapse are obtained as a special case. 
For a particular model we demonstrate that the radiating fluid sphere collapses 
without the appearance of the horizon at the boundary. 
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1 Introduction 

The study of radiating stars in the context of general relativity has generated 
much interest in researchers because of the variety of applications in relativistic 
astrophysics. These studies are important as they enable us to investigate physical 
features such as surface luminosity, dynamical stability, particle production at the 
stellar surface, relaxation effects, causal temperature gradients and other thermo- 
dynamical processes. Some relevant references investigating these issues are given 
by Di Prisco et al. pQ, Govender et al. [2J, Herrera et al. [3] and Pinheiro and Chan 
[2]. Relativistic radiating stars are also important in the process of gravitational 
collapse, describing the final state of stars, formation of singularities and black 
hole physics, in four and higher dimensions. Recent investigations in this regard 
are contained in the works of Goswami and Joshi [5 , Joshi [6] and Madhav et al. 
7 . In particular, the validity of the cosmic censorship conjecture can be tested in 
this physical scenario. 
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The model of a relativistic radiating star undergoing dissipation was completed 
by Santos [8] by analysing the junction conditions at the stellar surface. By match- 
ing a shear-free interior spacetime to the radiating Vaidya exterior spacetime, he 
showed that at the surface the pressure is nonvanishing and proportional to the 
heat flux. Subsequently several explicit relativistic radiating stellar models have 
been found by investigating the appropriate boundary condition. Kramer [9] and 
Maharaj and Govender [10] generated nonstatic radiating spheres from a static 
model by allowing certain parameters to become functions of time. Kolassis et al. 
11 and Thirukkanesh and Maharaj [12] assumed geodesic fluid trajectories to pro- 
duce new radiating models. In the approach of De Oliviera et al. [13] and Nogueira 
and Chan [14] the model has an initial static configuration before the radiating 
sphere starts gradually to collapse. Exact solutions for shear-free interiors which 
are conformally flat generate radiating stellar models as shown by Herrera et al. 
[15] , Herrera et al. [16] , Maharaj and Govender [T7] and Misthry et al. [18] . Stellar 
models which are radiating with nonzero shear are difficult to analyse because of 
the complexity of the boundary condition. However even in this case there have 
been advances in obtaining exact solutions. Particular exact models have been 
found by Naidu et al. [19] . Rajah and Maharaj [20] and Pinheiro and Chan [21] , 

In this paper we seek to generalise the Santos junction conditions by matching 
a shear-free interior spacetime to the generalised Vaidya exterior spacetime. The 
energy momentum tensor of the generalised Vaidya spacetime may be interpreted 
as a superposition of two fluids, a null dust and a null string fluid. The physical 
properties of the generalised Vaidya spacetime have been discussed by Husain [22] 
and Wang and Wu [23] . Glass and Krisch [2l] have interpreted the exterior space- 
time as a superposition of two fluids outside a relativistic star, the original Vaidya 
null fluid and a new null fluid composed of strings. By assuming diffusive trans- 
port for the string fluid Glass and Krisch [25] found new solutions to Einstein's 
equations with transverse stresses. Physically reasonable energy transport mech- 
anisms have been generated by Krisch and Glass [26] in the stellar interior with 
the generalised Vaidya metric as the exterior spacetime. These investigations, and 
other treatments, have largely focussed on physical processes in the exterior of the 
stellar model with a generalised Vaidya atmosphere. To fully describe a radiating 
stellar model requires generation of the junction conditions at the stellar surface. 

We follow the convention that the coupling constant and the speed of 
light c are unity; the metric has signature (— + ++)■ In Sect. 2 we present the 
field equations for the interior and exterior spacetimes. In Sect. 3 the matching 
of the interior and exterior spacetimes across the stellar surface is outlined. The 
new set of junction conditions are derived for the generalised Vaidya spacetime. 
We indicate how the new junction conditions generalise the junction conditions 
previously derived by Santos [8]. The physical significance of our new result is 
highlighted in terms of a string fluid. We also consider the new junction condition 
in the context of conservation of momentum flux across the stellar boundary. In 
Sect. 4 we generate a particular model and demonstrate that the new generalised 
junction condition has a solution. In this model it is possible for the radiating fluid 
sphere to collapse without the appearance of a boundary. The results of this paper 
are briefly summarised in Sect. 5. 
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2 Field Equations 

Spacetime needs to be divided into two distinct regions, the interior spacetime 
A4~ and the exterior spacetime _M + for a stellar model. The boundary of the star 
E serves as the matching surface for M~ and M + . The boundary or stellar sur- 
face is a timelike three-dimensional hypersurface. We require for the first junction 
condition that 

(ds + )s = (ds_)s = ds s (1) 

so that the line elements match on the boundary S. The second junction condition 
is generated by the continuity of the extrinsic curvature of S across the boundary 
given by 

(K+^)j: = (K- p ) s (2) 

Note that the junction conditions (fTJ) and ([2]) are equivalent to the Lichnerowicz 
[27] and O' Brien and Synge [28 junction conditions. 

The line element for the interior manifold .M - is given by 

ds 2 = —A 2 (t, r)dt 2 + B 2 (t, r)[dr 2 + r 2 {d0 2 + sin 2 9d<f> 2 )] (3) 

in comoving and isotropic coordinates. The interior spacetime is expanding and 
accelerating but is shear-free. A physically relevant interior matter distribution 
that is consistent with ([3]) is given by 

T ab = (A* + P) u aUb + pgab + q a U b + q b U a (4) 

where \i is the energy density, p is the isotropic pressure, q a is the radial heat 
flux vector and u a = ^Sq is the comoving fluid four-velocity. The Einstein field 
equations G~ b = T~ b for the interior manifold A4~ are given by 

B 2 1 ( B" B' 2 4B'\ 
M = 3^--^2— -_ + -_] Ca) 

1 ( „B B 2 „AB 



A 2 [ 2 B B i +2 AB 

1 I B' 2 n A' B' 2 A' 2B' \ 

B A B B 2 I A' 
2 — -7T— — + 2—- - — — 7 + 



A 2 B A 3 B A 2 B 2 r AB 2 

+ r B3 + AB 2 B A + S3 ( } 

2 ( B' B'B A' B\ 
q= -AB^\-^ + -B^+ AB) (5d) 

where dots and primes denote differentiation with respect to t and r respectively. 
The results (HJ-Q were first obtained by Santos [8]. 

The line element for the exterior manifold A4 + is taken to be 

ds 2 = - (l - 2 m( "' r) ) dv 2 - 2dvdr + r 2 {d8 2 + sin 2 Odcf) (6) 
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where m(v, r) is the mass function, and is related to the gravitational energy within 
a given radius r (Lake and Zannias [29], Poisson and Israel [30]). This metric is often 
called the generalised Vaidya spacetime since it reduces to the Vaidya spacetime 
when m = m(v) which is the mass of the star as measured by an observer at 
infinity. It has been demonstrated by Husain [22] and Wang and Wu [23] that an 
energy momentum tensor consistent with (|6]) is 

rp+ _ rp(n) rp{m) 

± ab — ± ab ' 1 ab 

= elj b + (p + P) (l a n b + l b n a ) + Pg ab (7) 

which represents a superposition of a null dust and a null string fluid. In general 
Tab represents a Type II fluid as defined by Hawking and Ellis [31] . The null vector 
l a is a double null eigenvector of the energy momentum tensor T^ b . The weak and 
strong energy conditions, and the dominant energy conditions are satisfied for 
proper choices of the mass function m(v, r). In ([7]) we have introduced the two null 
vectors 



7 rO _ 1 

ia — n a — ^ 



1-2^1 



Sa+Sl (8) 



where l a l a = n a n a = and l a n a = —1. The Einstein field equations G^ b = T^ b 
for the exterior manifold M + are then given by 

e = -2^ (9a) 



P = 2^ (9b) 



P = (9c) 



r 



where we have used the notation 



dm dm 
m v = — — , m r = — - 

av or 

We interpret e as the density of the null dust radiation; p and P are the null string 
density and null string pressure, respectively. Note that (T7))- ([9c)) were derived by 
Wang and Wu [23]. 



3 Generalised Santos conditions 

It is possible to match the spacetimes (|3]) and ([6]) across the boundary U. Since 
the derivation is similar to the Santos [8] treatment we provide only an outline of 
the argument for our more general case with m = m(v, r). The intrinsic metric to 
the hypersurface S is defined by 

ds% = -dr 2 + y 2 [dO 2 + sin 2 6d<f> 2 ^ (10) 
For the interior spacetime MT we obtain 



A(rz,t)dt = dr 
r s B(r s ,t) =y{r) 



(11a) 
(lib) 
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For the exterior region A4 + we generate the results 



i _2m +2 dr\ = /dv 
r dv J r \dr 



(12a) 
(12b) 



Equations (|11[) and (|12[) correspond to the first junction condition ([!]). Observe 
that the quantity r was defined on the surface £ as an intermediate variable. On 
eliminating r we have 



A{r E ,t)dt 

r s (v) = rB(r E ,t) 



, 2m „dr z \ 1/2 , 

1 h 2—^ ] dv 

r-£ dv 



(13a) 
(13b) 



Equations (|13[l are the necessary and sufficient conditions for the first junction 
condition ([TJ) to be valid. 

The intrinsic curvature for the interior spacetime A4~ has the form 



AT 



AT 



B A 
~{rB)'} 



A 33 = sin 9K 22 

The extrinsic curvature for the exterior spacetime A4 + has the form 



Kt 



V _m _TO r 
- - V^r + V 

v r 



v 1 



2m 



r + rr 



A"+ = sin 2 0A^ 



(14a) 

(14b) 
(14c) 

(15a) 

(15b) 
(15c) 



where r = ^ and l> = 4r. Observe the appearance of the term containing m r 
in A^ which does not exist in the treatment of Santos [8]. As we shall see later 
this has a profound effect on the physics of the model. Equations (|14[) and (115|) 
correspond to the second junction condition ([2]). The mass profile in terms of the 
metric functions can be generated by eliminating r, r and v. We observe that 



i(t),r) 



rB 
~2 



1 + r z 



,B^_ 

A 2 



(16) 



which is the total gravitational energy contained within the stellar surface S. We 
also establish the relationship 
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This expression may be simplified further: multiply with 1 + r-g- + r-j and utilise 
(|5b)l and (I5rf|) . We then arrive at the result 



P={«B-^) s (18) 



which generalises the junction condition of Santos [8]. Hence we have demonstrated 
that the junction conditions ^ are equivalent to 



Equations (|19p are the necessary and sufficient conditions for the second junc- 
tion condition ([2]) to be valid. We point out that the mathematical approach and 
procedure that we have followed is similar to Santos [8] . However in our case the 
external stellar atmosphere is the generalised Vaidya spacetime. The form of the 
equations (|10p - (|19b|) can be related to the equations of Santos [8] since they have a 
similar structure. The fact that m = m(v, r) fundamentally affects the final result. 
The equations for the extrinsic curvature ()15of) . the matching condition (|17[) and 
junction condition (|18[) are fundamentally different. 

We have generated the relationships (|13|) and (|19[) so that the junction con- 
ditions ([1} and ([2]) are satisfied for the shear-free interior spacetime ([3]) and the 
generalised Vaidya exterior spacetime ([6]) across the hypersurface S. This gener- 
alises the Santos [8] result for a relativistic radiating star when m = m(v). Observe 
that when m depends on the coordinate v only then (|19b|) becomes 



at the boundary S, which is the earlier Santos junction condition. When (|20[) is 
valid then the pressure p on the boundary depends only on the heat flux q. We 
have shown here that if m = m(v, r) then (|19fo[) is valid, and the pressure p on the 
boundary depends on the heat flux q and the gradient m r (v, r). 

The generalised Vaidya spacetime has physical significance and contains many 
known solutions of the Einstein field equations with spherical symmetry. It con- 
tains the monopole solution, the de Sitter and Anti-de Sitter solutions, the charged 
Vaidya solution and the radiating dyon solution. The physical features and the en- 
ergy momentum complexes, that provide acceptable energy momentum distribu- 
tions for these systems, have been studied by Barriola and Vilenkin [32], Chamorro 
and Virbhadra [33], Virbhadra GI]-[36] and Yang [37]. Glass and Krisch [23]-[25] 
and Krisch and Glass [26] have interpreted the generalised Vaidya spacetime to 
represent a superposition of an atmosphere composed of two fluids: a string fluid 
and a null dust fluid. This atmosphere may model several physical situations at 
different distance scales, eg. the exterior regions of black holes (distance scale of 
multiples of the Schwarzschild radius) and globular clusters containing a com- 
ponent of dark matter (distance scale of the order of parsecs). The additional 
term 2^ in the boundary condition (|196[) arises from the matching at the surface 
S. This quantity has physical significance and can be interpreted as a particular 




(19b) 



(19a) 



p = qB 



(20) 
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contribution from the energy momentum tensor. We observe that the term 2^ in 
(j 19£>|) is the same quantity as that in . Therefore we may interpret the quantity 
2^ as the string density p. 

We can therefore write (|19b[) in the more transparent form 



at the boundary S. Consequently for a radiating star with outgoing dissipation in 
the form of radial heat flow, with the generalised Vaidya spacetime as the exterior, 
the pressure on the surface depends on the interior heat flux q and the exterior 
string density p. The appearance of the quantity p in (|21[) allows for more general 
behaviour that was the case in the Santos [8j treatment. From (|20[) we observe 
that q = implies that p = on £ and the exterior manifold M + must be the 
Schwarzschild exterior metric with m being constant. In (121[) we note that we 
obtain the Schwarzschild exterior geometry when q = = p. However it is clear 
from (|2ip that when qB = p then p = on S and the exterior spacetime remains 
the generalised Vaidya spacetime with m = m(v, r). In addition, when q = then 
p = —p on S and the interior is not radiating. 

It is possible to provide a physical interpretation of our result by consideration 
of the momentum flux across the boundary S. Since the quantity ()19a[) represents 
the total gravitational energy for a sphere of radius r within £ we can write 
m(v , r) = m(t, r). Partially differentiating (|19al) eventually leads to the result 



which reduces to the corresponding Santos [8] equation when m = m(v). The 
radial flux of momentum across the hypersurface £ is defined by 



where e a and n are vectors which are respectively tangent and normal to S. 
For conservation of momentum flux across £ we must have 



p=[qB- p] s 



(21) 




(22) 



e n T ab 



(23) 



F + 



= F~ 



(24) 



In the interior manifold M we can generate the quantity 



F~ = — 



qB 



(25) 



In the exterior manifold A4 we can produce the quantity 




2 -2 

— V 777, 



"V 



(26) 



Then equations ([22j) - (J26J) yield the result 




which is the same as (|19b|) . Therefore the junction condition (|196[) corresponds 
to the conservation of the radial flux of momentum across the hypersurface S. It 
represents the local conservation of momentum. 
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4 A particular model 



To illustrate the utility of the generalised Santos condition we consider a specific 
example relating to horizons. For simplicity we take the mass function to be of 
the form m = pr 2 where p is an arbitrary constant. For this example we consider 
a particular form of the metric coefficients given in . We choose the coefficients 
to be separable in r and t so that 



A = a(r), B = b(r)R(t) 

The field equations ([5]) then yield the quantities 



1 

R? 

J_ 

R? 



—^R 2 



1 2b" 



b 2 



a 

2a' R 



1_ / 

~2 



2RR + R 2 ) + 



b^_ 4V_ 

b 2 rb 

1_ 

62 I b 2 



+ 



2a' V 



ab 



+ 



a o 

- + T 

a o 



R 3 a 2 b 2 

and the condition for the isotropy of pressure 



(,'2 



-J If 1 . 

^ + b --2 b - 

a b b 2 



I it 

} a o 
ab 



ra rb 



The boundary condition (|19&|l now yields at r = Te the equation 

2RR + R 2 + MR = M + V 

where we have used the particular choice of the mass function m 
potentials (|27[) and the system (|28p . In the above we have set 

ka 2a ' 
M = -- 



Af - 
V 



a 
b 2 



b 2 ab ^ r 



a o 

- + T 

a o 



, a 2 p 
' b 2 



(27) 

(28a) 
(28b) 
(28c) 

(29) 

(30) 
pr 2 , the 

(31a) 
(31b) 

(31c) 



which are constants in the integration process as (|30p holds on the boundary S. 

Equation (|30[) does not have a general solution in closed form; for the purpose 
of our investigation we observe that it admits an elementary particular solution of 
the form 



R(t) 



-Ct 



(32) 



where C is given by 



M± [M +4 
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This solution is useful in that at the surface of the collapsing star the ratio 2L£ 
independent of time. The mass profile given by (|19ap and the solution (|32 
the following expression 



is 

yield 



2mi 



2 b 2 C 2 
2a 2 



rb 
T 



r 2 b'- 



2b 2 



(33) 



which is valid on the stellar surface. It is clear that the quantities in the above 
equation are evaluated at the boundary E; there is sufficient freedom provided by 
r, a, b and C so that the right hand side can be constrained in such a manner that 



the ratio 



is strictly less than unity. Consequently no horizon will form. We 



can illustrate this explicitly by setting b(r) 
£ are arbitrary constants. Then we get 



2m 5 



1 and a(r) = \^r 2 + (3 where f3 and 



C 2 r 2 



(34) 



We choose C 2 to be less than ^ +/3£-r4£ 2 r 2 on U so that 1- 



is not zero for all 



2m z 

4 ^ / un ^ nu uiiau ± — 

time. Consequently the boundary surface never reaches the horizon. This simple 
example demonstrates the absence of the horizon and is similar to the result of 
Banerjee et al. [38] which holds in the conventional Vaidya spacetime. We regain 
their results when V = in (|3ip . We interpret our result to mean that there is no 
accumulation of energy in the interior of the radiating star as it releases energy at 
the rate of generation during the collapse. 



5 Conclusion 

In this work we have produced a general model of a relativistic radiating star by 
performing the smooth matching of a shear-free interior spacetime to the gener- 
alised Vaidya exterior spacetime, across a timelike spatial hypersurface. We have 
demonstrated that with the generalised Vaidya radiating metric, the junction con- 
ditions on the stellar surface change substantially, and consequently represents a 
more general atmosphere surrounding the star. The atmosphere is a superposition 
of null dust and a string fluid. We find that the density of the string fluid affects 
the pressure at the stellar boundary. We have shown explicitly that 

V = QB — pstring 

at the stellar surface. If the weak and strong energy conditions or the dominant 
energy conditions are satisfied then pstring > (p / 0) and p s tring > Pstring > 
(p ^ 0) respectively. This indicates that for outgoing heat flux in gravitational 
collapse, the exterior string fluid density reduces the pressure on the stellar bound- 
ary. It is interesting to note that we have shown using a geometric approach that 
the derivative of the mass function with respect to the exterior radial coordinate 
is related to the string density. By means of a simple example we have considered 
the formation of horizons for a heat conducting sphere which radiates energy into 
the generalised Vaidya external atmosphere during collapse. The horizon does not 
appear in this example at any stage of the collapse. 
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